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During the period of July, August and September, 1981, progress has 
continued on all phases of the research program. Particular attention has 
been given during this period to the preparation of papers and documents 
related to the Magsat program. 

On July 30 and 31, 1981, we participated in the Workshop for Magsat 
Investigators held in Edinburgh, Scotland. Following these workshops, we 
presented two papers on the 4th of August 1981 in Session II entitled "Scien- 
tific Results from Magsat" at the Fourth International Association of Geo- 
magnetism and Aeronomy Scientific Assembly. The abstracts for these papers 
are published in lAGA program for the Fourth Scientific Assembly. The first 
paper presented, entitled “US Aeromagnetic and Satellite Magnetic Comparisons" 
by R.R.B. von Frese, W.J. Hinze, J.L. Sexton and L.W. Braile, made a detailed 
comparison between the upward continued NOO aeromagnetic and POGO magnetic 
anomaly data over the conterminous U.S. The second paper, entitled "Satellite 
Elevation Magnetic and Gravity Models of Major South American Plate Tectonic 
Features" authored by M.B. Longacre, W.J. Hinze, R.R.B. von Frese, L.W. 

Braile, E.G. Lidiak and G.R. Keller presented the results of our recent analysis 
and processing of the Magsat data over South America, the Caribbean and adjacent 
marine areas and the compilation and analysis of the associated gravity and 
other geophysical and geological data. The results are enumerated in the 
abstract and include modeling of the Andean subduction zone, the continent 
of South America and the Amazon River Aulacogen. 

We have also prepared several figures and related discussions for documents 
being prepared by NASA on the feasibility and applications of the Gravsat 
Mission and the proposed Magsat B Mission. In addition, we have completed 
the preparation and review of a manuscript on "Long-Wavelength Aeromagnetic 
Anomaly Map of the Conterminous U.S." which is in the process of being sub- 
mitted for publication. 

We have also received confirmation of the publication of a paper by 
R.R.B. von Frese, W.J. Hinze and L.W. Braile on "Long-Wavelength Anomaly 
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Map of the U.S." in the Geophysical Journal. In addition, the paper by R.R.B. 
von Frese, W.J. Hinze, L.W. Braile and A.J. Luca entitled “Spherical-Earth 
Gravity and Magnetic Anomaly Mod'*ling by Gauss-Legendre Quadrature Integration" 
has been published by the Journal of Geophysics (a copy of this paper is 
enclosed). 

Current activities are centered around the preparation of papers for 
the forthcoming special issue of the Journal of Geophysical Research Letters. 
Papers are being prepared related to our continuing studies and analysis 
of the U.S. Naval Oceanographic Office Magnetic Survey Program of the U.S. 
and our Magsat studies in South America. 

Total expenditures as of October 1, 1981 are approximately $58,511. 
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Sphcrical-tlarth Gravity and Magnetic Anomaly Modeling 
by Gauss-I^endre Quadrature Integration 

R.R.H. vtiii I'rcsc'. V\..l. llinzc', 1 W. I)r;iilc‘. ;iiu! A.J. l.iica’ 

' IX-[virimci)i of (ioiMiciu.vs. 1‘tinluc I'liixorsii). Wtsi 1 ura>citc. Iiuliiimi 47907. l^SA 
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Absirtiri. I lie iiv.-riilncsN of toii(;-wavdci>i!lh p<’tcntial field 
un«>inalies m liihospherii.' iiiierprelation is pieatly inercased 
with splicrieal earth niodeliiig leehiiiqucs. ( iaiiss-i.egeiidrc 
quadrature integration used to calculate the anomalous po- 
tential *>f gravity and magnetic fields and their spatial de- 
rixativcN on a spherieal earth for an arbitrary boily represented 
by an equivalent point source distribution of graxiiy (Soles or 
magnetic dipoles. I he distribution of equivalent point sources 
is de'crniined directly from the coordinate limits of the source 
volume. Variable integration limits for an arbitrarily shaped 
biHl) are ilenved from interpolation *>f pxiinis which approxi- 
mate the body's surface envelo|x*. Ihe versatility of the method 
is enhanced by the ability to treat physical pro|X‘ity variations 
within the source volume and to consider variable magnetic 
fields over the somce ami observation surface. A number of 
examples verify and illusiiate the capab.lities of the technique, 
including preliintiiary modeling of potential field signature's for 
Mississippi embaymeni crustal structure at satellite elevations. 

Key wwds: Spherical coordinate gravity ihodeling -Spherical 
coordinate magnclic iixHleling I quivalcnt sources -(iatiss- 
l.cgendrc qiiadraltire miegralion Satellite jiotetitial fields 
- Mississippi embaymeni 


Introduelion 

C'lvmpulation of theoretical anomalous gravity and magnetic 
fields from geologic models is an impvirlaiU element in in- 
terpreting potential field data and designing surveys, l or re- 
gional giavity and magnetic surveys measured in degrees v>f 
latitude and longitude, procedures are desired which model 
directly, in spherical eoordinaies. potential field anomalies due 
Itv large-scale, arbilrarily shaped sources of variable density or 
magiieti/alivtii chaiaelerislies. 

In a review of compuier modeling techniques Hhaltacluiryya 
(I97X( presenis meliiods that in principle are suitable for 
spherical earth modeling of regional features. I hese procedures 
are generally based on approximations of ihe anomalous 
source as a group of prisms or polygonal laminae, the elTeets of 
which are evaliialed and summed at each observation point to 
yield the total anomaly However, for typical spherical earth 
m» ‘ ' 'g applications, the book-keeping problem involved witli 
siihvtividiiig the laige-se.ilc soime into simple loiiiis to rellcct 
aihiliaiy cliaiacleiishis <•! gemueliy ami physic-d pio|K'iiies is 
eonimonly lormidable 

A simple and more cITicicni procciluie is to lepreseiil 


quadratures of Ihe source volume by equivalent p«>int sources 
according to the well known leehniqiie of (iattssian quadrature 
integration. In pri iciple. the appropriate geometrie distribution 
of cjuivalcnt point source's can be determined directly from the 
coordinate linrils of the source volume, so that an accurate 
estimate *'f the source-afliliatc'd anomaly is obtained by 
evaluating and summing at each observation point the anomaly 
values due to each point source of the equivalent source distri- 
bution. 

Ciaussian quadrature is a timo-bonored teehnique for 
numerical integration and has been well studied in the litera- 
ture of numerical methods (Carnahan et al. I%‘>l. rolcnlial 
field nuHleling by (iauss-l.cgemlre quadrature was used by Ku 
il‘)77) to evaluate gravity and magnetic anomalies in Cartesian 
coordinates due to bodies of arbitrary shape aiul magnetic 
pvilari Zillion. In this discussion, the method is extended to 
spherical coordinates and the general problem of sources with 
arbitrary sbape and variable density and magnetization proper- 
ties. 


Description of Method 

As illustrated in I'ig. I. it is eonvenient for gravity and mag- 
netic nuHleling problems to c'onsider the anomalous body as 
being composed of a source volume distribution of gravity 
|Himi poles (ir magiietie point dipoles, respeelively. lienee, lo 
estimate the anomalous gravity or magnetic field at some ob- 
servaiion point, it is nccx'ssary lo evaluate and sum the ano- 
maly values due to each of these point sources at Ihe obser- 
vation point. 

In particular, the radial anomaly. <tg. due to a gravity 
point pole referred lo a geoecnlrie coordinate system as shown 
in I ig I is given by 



where 


a - universal gravitational constant (=6.h7xl0 " enr'/ 
(.x'se'e*)), 

ft =diMaiice from observation point {r.0.if>) In source point 
{r'-O'.ifi'l 

i.r' ladial distances from earth's eenler lo the observalioii 
l>t»iiil and source isoiiil. repeeitvely, 

II. I) eii-laliliide eiHtrdinales of observation and soiirtx* points. 
res|vclivc'ly. 
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i/i. I'l' - IoiiuiIikIv i‘Oi>i\liii;iles nf oh\i-rv;ilioii ;iiul M'ura’ poiiils. 
ii'\|X'cti\cl\. ;iiul 

till - iiuiNS coni rust of ilic iioiiii pole, 

lor n magiiL'iic point dipole, on tlic olhet luiml. ihc total 
nuitiK’iic iniciiMiv aiionial). . 1 1. ii> jiiMn h> 

'' 1 “ ‘ (“ ‘"(k))} ■' 

where 

I, r -= gradient o(X'ialt)r,s in ob'<ei\alion point and Miniee point 
Cl>ordlnale^, res|X’eli\el> 

1/ - magneii/aiiitn eonira>i of the |ioini dipole, and 

u u iinil geomagnetic field \ecliM> at the olKei\aiion point 
aiul soniee itoinl. lespeciivel) 

( omeniionalU. the nnil vecloi> u and u are e\pre>sed in terms 
of geoinagnelic field nicliiuiion (/./ ) and declination {1)1)1 
\lso. when the pt'ini dipole |solari/atioii is h> mdiiclion 

It lA / (.'1 

w here 

lA - niagnelic suscepiibilil) coniiasi of the point dipole, aiul 
I scalai geoni.igneiic lielil iiilensiiv .it the souict' ponii 

Hence, in regional-scale applic.itioiis ol I i.|. (.1 geoinagnelic 
lield models '.ncli as the UiKI l‘*(>'' (( am el .d I‘hi7| ,ne 
noiiiiall\ iiseii to obtain peilineni xalnes ol (/./>./ I at the 
source point and (/. I){ at the observation point 


Consideration of P.qs. (If and |2| shows that the incremental 
gravity and tiiagnclic anomaly values due to a point source 
located a distance R from the observation point can be general- 
ized according to the relation 

</(/<) /»v (41 

where 

!/(/() = the geometrical point source runction which describes 
the inverse distance between the ohscrvalion point and 
source point, and 

d.v =the appropriate physical property contrast of the point 
sourcs:. 

Hence, to determine the total potential Held anomaly it is 
necessary to evaluate at each observation point (r. 0. </>) the 
volume integral given by 

Wt. r'l, 

j J j (</(r, (A</’;r.f). r/)').4.v). illl'tltji' (.‘'I 

where the primed variables refer to the coordinates of the 
anomalous body such that 

i/r. (/i|, = lower and upper longitude limits of the source volume. 

= lower and upper co-latitude limits of the source vol- 
ume. and 

r„.r'„ = lower and upper radial limits of the source volume. 

Consider, now. the numerical evaluation of the innermost 
integral in liq. (5l. Most numerical integration techniques in- 
volve the use of interpolation polynomials to approsimate the 
integrand according to a summation formula of the general 
type given by 

r». Hi 

I (</(r’M \)(/r‘ i V ,'t,<i(r,) .1 V (f>( 

I I 

where the iii values .4, arc the weights to be given to the ui 
functional values (</(rj) J.v) evaluated at the interpolation cinir- 
dinutes, r'. Conventional integration formulae such as 
Simpson's rule or the trapezoidal rule assume equal spacing of 
the arguments r' which is generally appropriate when dealing 
with an integrand that is not well known analytically. However, 
as shown in Iqs. (I) and (2| the integrand being eonsidcrcvl 
here involves a familar analytic function which may be com- 
puted for any argument to great precision. In such instances, 
(i.mssian quadrature formulae can he developed to yield selec- 
ted values of interpolation r| and coefricients 4, so that the sum 
in liq. (ft) gives the integral exactly when (</(r ).1x) is a pi*ly- 
nomial of degree 2»i or less (Carnahan el al. 1%')). 

In general, it may be shown- that the (iimssian cocITieieuts 
I, can be obtained from a polynomial of order iii which is 
orthogonal over the interval of integration such that the ni 
(voints of interpolation r| are the zeros of the polynomial. 

I amilies of orthogonal polynomials which are commonly used 
to develop (iaussian quadrature formulae include l.egcndre. 

I ague I re. Chebyshev and Hermite pvilyiiomials. However, in 
this discussion, only the prototype of the Gaussian method in- 
volving I.egendre polynomials is considcretl. 

I egendie |>olynomials /’,(f'l of order m which are ortho- 
gonal over the interval - I g I are given by 

(2-L!)Cr"-’^"’‘ 
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Ac«:ori]in):ly. the slandartl (iatiss-l.cgciulrc quatliaUirc over ihc 
interval ( 1 . 1 ) is giu-ii li> 

I m 


f («/(r).lv|i/r‘ ' Y. (S) 

I I I 


where the interpolation points F, at whieh the integraiui is 
evaluated are the /ems of r\|. (7) and the Ciaiissian eoeffieienls 
arc 




2 ( 1 -r-;-) 




Now for arbitrary limits of integration sueh as indicated it) 
Lq. ( 6 ) it is necessary t») map the standard interval - 1 <r;g 1 
into the interval of integration r],gr|grj, according to the 
transforniation 


• ■) 


( 10 ) 


Accordingly, the integral in I q. ( 6 ) can be appioxinialed us 


Y (,</(r;).iv 
^ 1.1 


(li) 


hxtending this procedure to volume integrals is straightfor- 
ward. Thus, the ('lauss-l.egendrc formula for the general eva- 
luation of I q (.S) is given by 




I j j O’, <j)'i /f sjdr'Jlf' Jif)' 
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where 

0.5 |r;tr;, r;“) ( t r;..| 

^1=0.5 {<?;(«/>;»- 1/»;,) + </>u + </>;j 

= coordinates of the subdivision in the limits 'f in- 
tegration from - I to I which e*>rresiroi)d to zero 
nodes of llq. (7). 

A,. /lj = (jauss-Lcgendre quadrature coefficients given by I'q. 
(9). 

-lower and upper longitude limits of the body for the 
A-th longitude eomp*»nent of tlie equivalent point 
source coordinates. 

-lower atid up|vr eo-latitude litnils *if the body for 
(he j-th eo-laiitude compvment i>f the equivalent 
point SI lice eooidmates. and 
lower and uppet ladial limits of the body for the 
i-th radial component of the equivalent point souiee 
coordinates. 


The quadrature formula given in I q (12) shows that gravity 
and magnetic anomalies can lie computed accurately by sum- 
ming at each observation point the anomalous efTeet of nk > m 
X III ei|uivalent |ioml souivvs li>cale«l at souiee iwiiiit coot 
dinates (r|.()]. 04 ). where each of the differential point source 
anomalies is appropriately weighted by (iauss I egendre 


quadrature ciK'lfieients and the volume coordinate limits of the 
anomalous body being nuHleled. 1 'his result is well suited to 
machine computiilion where the input consists principally of 
values of the iuiegrand for selected sivurcc points (r], (l^.iA^I. 
affiliated subdivision ciHtrdinutes (C'. ^i), and coclficients 

/lc)< ttnri the volume coordinate limits of the body for 
each dimension of every source point crnvrdinate (r'.fl'^.^i). 

1 he selected values of the integrand for gravity or magitetic 
modeling purposes are readily obtainc*d from liqs. (I) or (2). 
respectively. Also, (he Legendre subdivision coordinates of the 
interval ( I. I) and associated Ciatissiun coclficients may Iv 
computed directly from l-qs. (7) and (9). respectively. However, 
it is generally found to be more machine efneient to input these 
values from tables using an algorithm such as dcsc'rihcd hy 
Carnahan e( ul. (I969). Values which arc applicable for such 
algorithms are tuhululed to .fO digit precision for orders n = 2 
5I2 by Stroud and Seercsl (1966). However, for mos« regional 
lithospheric modeling applications expcricnc'c suggests that sub- 
division coordinates and associated coefficients to lO digit pre- 
cision for orders up to ri = I 6 arc normally sufTidcnt. 

For a uniformly dimensioned body such us a prism, the 
integration limits for the evaluation of Hq. ( 1 2 ) arc easy to 
sjK’cify. In this case, for example. (0l„. l/>^^) = (</'l. («„,. 

=(((„. (^(,) and (r'„,r'^) = (r;, .*•;). However, for the more general 
case of a body with arbitrary shape the integration limits arc 
known in ooe dimension only, so (hat the problem of determin- 
ing the integration limits in the remaining two dimensions for 
each equivalent point source coordinate must be eoiisklcrcd. 

Procedures cun he developed to handle this prvtbleni ef- 
ficiently, such as Ihc method described by Ku (197?) where a 
modified cubic spline funclion is used to interpolate the desired 
integration limits from a set of body point coordinates which 
provide a rough approximation of Ihc surface envelope of the 
body. Typically, the procedure is to specify, for example, the 
longitudinal limits of integration of the hmly to obtain the uk 
(iauss-l, egendre iuhIcs (/>» as described above. Interpolations ol 
the lw>dy p*)inl coordinates arc performed next to determine 
the maximum and tninitnum latitude coordinates of the body 
for each longitude cvHirdinalc ^ 4 . These values t'f e<>urse pro- 
vide Ihc latitude limits of integration for evaluating the m 
nodes (f^. Similarly, the radial coordinates of the body points 
arc interpolated ul each horizontal coordinate (<!>[. f)^) to yield 
appropriate radial limits of integration from which Ihc til nodes 
r' can be determined. Procedures such as this arc readily adap- 
ted for efficient machine processing so that the integration 
limits of arbitrarily shaped bodies can be determined accurately 
for evaluation of the quadrature formula given in Eq. (12). 

The quadrature formula in general has considerable versa- 
tility in modeling applications tecause anomalous gravity and 
niagnetic potentials and their respective sptitiul derivatives t>f 
any order are all linearly related. Hence, to model the r.adiul 
derivative of the potential field anomaly due to an arbitrary 
source, for example, it is necessary simply to exchange <i{R) for 
in liq. (12). Additional geophysically interesting quan- 
lilies which can he modeled from simple linear trumsforntalions 
of the integrand of Eq. (12) include the anomalous potential, 
vector anomaly components and the spatial derivatives of any 
order. Relative go'idal anomalies can also be motlclcd hy 
computing the anomalous gravitational potential of the body at 
the stirface of Ihe earth and dividing il hy normal gravity 
.iccoidiiig to lliun's Ibiiiiiila (llciskanen and Mori)/ l%7). 

Itudics with variable physical property contrast arc also 
accomm«Hiatcd readily by the quadrature .solution. To cmpha.s- 
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FIk. 2. Ciruviiy unomuly 
Gompari!ion!> for a spherical prism 


izc ihc fad lhal /l.v can be c.xprc.s,scd a fundion of source 
position (i.c., = (#■.(/>')) when necessary, the gcncrali/cd 

point source anomaly (i/(K).^1.\| has been carried intact 
throughout the foregoing developments. Conventional model- 
ing procedures normally assume bodies with uniform physical 
properties so that each property variation must be modeled as 
a separate body. To model the total intensity magnetic anom- 
aly due to a regioiuil source subjcd to regional variations of 
geomagnetic field induetion by (he well known method of 
I'alwani il9liM. for csampic. requires that the source be sub- 
divided into blocks wherein the geomagnetic Held polarization 
is uniform. l£ach bKxk in turn must be subdivided into a group 
of approximating polygonal laminae to ensure accuracy of the 
numerical integration which evaluates the anomaly. Uy con- 
trast. the Gauss- Legendre quadrature approach is simply to 
polarize each point dipole in the quadrature formula according 
to the polarization characteristics of an acceptable numerical 
model of the geomagnetic Held (e.g.. the IGRI -I'f65 by Cain 
et al. I%7). Ihis example alst* illustrates the technique for 
modeling IhhIics with remanent magnetization by quadiatuic 
formulation. Here, the magnetization of the point dipoles is 
achieved using a polarization field that represents the vector 
sum of induced and remanent magnetic polarizations. 

Finally, considerable computational flexibility is available 
for the practical implementation of the quadrature formulation 
Ku (1977) noted the trudc-olTs which incur in applications of 
(he method between efTicient computation and (he accuracy of 
the solution. A large numlK'i of Ciauss-I cgciidic nodes ensures 
a very accurate quadrature solution, although the same degree 
of accuracy cun often be achieved by a substantially smaller 
numher of eqimalent point sources. In fact, the accuracy of the 
solution remains essentially unchanged for difTcrent numbers of 
nodes as long as the niHlc spacing is smaller than (he distance 
to (he observation point. Hence, the minimum luimiKr of no- 
des specified in a given application should be such (hat tlw 
distance to the observation point is greater than the nude 
spacing. In getwral, then, the accuracy of the quadrature for- 
mulation can be readily controlled by adjusting the inimK'r of 
Cjauss-Lcgendrc nodes relative to the elevation of die obser- 
vation point. 


llliKiratioa of Mclliod 

A computer program is described in von Frese et al. (I'JKOI 
which was developed for regional lithospheric gravity and mag- 
netic modeling applications by Gauss- Legendre quadrature. 
The program was used to construct examples that illustrate 
some of (he capabilities of the method. 

To demonstrate and verify (he mcthixl for regional-scale 
modeling applications, gravity and magnetic anomalic's due to 
a three-dimensional spherical prism were ’"'Hided on a spatial 
scale small enough (hat (he results could oc compared readily 
with conventional modeling techniques in Cartesian coor- 
dinates. In Fig. 2, for example, a comparison is made between 
(he (wo coordinate systems fur calculations of gravity anomalies 
due to a 6.67 km thick prism with density contrast 0.25 g/cm^. The 
observation grid consists of (16,16) stations uniformly spanning 
the region (27.t 274) *E, (40 41 ) *N at an elevation (Z) of 6.67 km 
above the top of the prism. 

fhe well known method of Talwani and Fwing (F>fi0) was 
used to calculate the gravity effect of the prismatic model as 
shown in Fig. 2A. The gravity anomaly was determined by 
evaluating and summing at each observation point the gravity 
anomalies due to II horizontal polygonal laminae used to 
approximate the prism. As shown in Fig. 2A, the resultant 
anomaly has an amplitude range (AR) between 36.0 mgal and 
0..3 mgul and an amplitude mean (AM) of 5.4 mgal. The anomaly 
was computed air! contoured in ('artesian ciMirdinates assum- 
ing I " 100 km. Ihis assumption distorts slightly the (rue geo- 

metry of a spherical prism in this 1° x |° region in northwes- 
tern Indiana where I* in longitude or latitude is more nearly 
equal to KK km or 1 12km. respectively. 

In Fig. 2B the gravity anomaly of (he spherical prism com- 
puted by Gauss'Lcgendrc quadr.aturc integration is contoured 
on a stcrcographic equal-area polar (SliAP) projection. In (his 
case, Ihc gravity anomaly was calculated in spherical axtr- 
dinates by evaluating al each observation point an nk xnj xnf 
>8x8x2»l28 point quadrature formula where the inte- 
gration limits were specified directly from the spherical coor- 
dinate limits of the prism volume. The results shown in Fig. 2 
indicate that the (wo methixls agree very well with respect to 
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Ki|t..L MiiiiiKiic ;inimuil\ 

lor ;» pi'isin 


llic rdalivdy MOall Npalial M.aL- under (.on>idciaiioii. Ihc sli|!hl 
disparity hciwccii (ho results is probably duo it) ycomotrio 
distortion inlrodtiood by tho use of C'artosian ooordinatos for 
tailoulaiion of tho sphorioal prism anomaly. 

A ooinparisttn I'f inagnotio anomaly fields is illusirulcd in 
I'ig. 3 for (ho prism iuikIoI ‘AUh voliimo niagnotio susooptibilily 
contrast oqiial to (),Ot)(I.Soimi om * (4itx».OOO.S in .SI units) and 
uniform |Mi|ari/ing intoi'sity i WI.OOO gamin. i (<><), UIO n 1 1 . in- 
clination r = iy and iloolinalion I) -O''. Ovoi tho obsorsalion 
grid, imiform goomatinoiio fidd atliludo oliataoloristios of in- 
clination / = 75 .iml dodinaiioii I)- O' arc also assuniotl. I bo 
total intoio'iy niagnolio anoinaly of (bo spboiioal prisin oaloti- 
luted III ( artesian oiHirdinaios according to the molbi>d of 
Talwani |l%.3) is sbovui in I ig. .3A. I bo lostibani anomaly 
oorros|vtiids moII uiiIi (bo magnotio anomaly oaloulatod m 
spborioal ooordinatos by (iauss-l ogoiidro qtiadraiiiro in l ig. .Ml. 

As a mote rogional-soalo osainplo, oonsidor the application 
of (iaiiss-l ogentiro tiiiailratiiro integration for moiloling gr.ivity 
and magnotio anom.ily signatures, at the .satellite dotation ol 
4M)kin. duo to (be crustal tbicknoss anomaly ttbicb is griildcd 
in l ig. 4.A. Ibis /one of onbancod crustal thickness, o.xtoiubng 
from the Texas panbandio n« riboastward into Kansas, is por- 
trayed by scisinic ovirlcnoo iWarron and lloaly b>7.3) as roughly 
10 km of lighter crustal material displacing denser (mantle) 
material l») a depth of .<0km. l o illustrate the gravity modeling 
procedure, a density coiiliast of ()..3gcm' teas assumed for 
this feature. 

Magnetic eonsitleraiions. on the other band, suggest that 
the crustal thickness anomaly may represent a /one of positive 
magneti/ation contrast due to dossnssard deflection ol the ( u- 
rie isotherm in the region of thick cooler crust.d material or 
thicker m.ignetic erusi extending inio non-magnetic mantle 
material. I'or the puriisises of this example, a volume suscepti- 
bility contrast of 10 ’emu cm' vt.is assumed for the /one of 
enhanced crustal thickness vdiicli is representative of the genci- 
ul magneti/ation ic|>oried lor (he lower crust (Hall 1974; Shucy 
ct al. 1 97.3). 

To determine the variable limits of integraiion loi the eval- 
uatHin of (he qu.idtattiie loiimila. (he crustal thickness anom- 
aly was referred to the body point grid shown in I ig. 4A I vir 


each (iauss-l. egendre node in the </»-longitude coordinates ol 
the binly. the O'-lalilude limits of the b*)dy were interpolaterl 
from the 6 points describing the body's boundary. The radial 
limits, in turn, were interpolated from the 6 boundary points 
and .3 interior points of the body for each Gauss-l.egendre 
node with horizontal (</>', (>')-e»H>rdinatc’s. For the particular nuHl- 
el considered here, of eoursc. the radial coordinates of the 9 binly 
IHvinls Used to approximate the subsurl'acv conrignration of the 
bixly were speeilied to represent a unirorm thickness of III kin 

Hie resultant gravity anomaly al 4J0km elevation due to 
(his feature is illustraic'd in I'ig. 413. The gravity effeei was 
caleulaled by evaluating an nk x nj x iii = 16 x 16 x 2 ^ .^1 2 |xnnt 
(iauss-l egendre quadrature formula over ihc (41. 41) grid of 
observations. The magnetic effect tvf the crustal thickness anom- 
,ily was al.so computed in this manner at 450 km elev.ition 
The resultant total intensity magnetic anomaly in the KiKi - 
|9(is u|xla(ed to I9(>K is demonstrated m I ig. 5A. lo remove 
amnnaiy distortion due to regional variability of the reference 
lield, the magnetic elfeet was next eomputcxi. assuming a tiinfoiin 
polarizing field intensity of (lO.(KN) gamma and radial geomag- 
netic field inclination at Imth source and observation |vvinis 
Accordingly. I ig 513 illustrates the resultant magnetic aimmaly 
field reduced to radial polariziition at 450 kin elevation for the 
crustal thickness nuHlel. 

l inally, for accurate implementation of the method the 
distance between the cx|uivalent sourex* points and (he obser- 
vation (xvini must be greater than the ev(uivalen( sourc-c point 
spacing within (be body. This limitation can be minimized in 
practice by either subdividing the body into an appropriate 
number of smaller bodies, or increasing the number of equiva- 
lent point sources, or increasing the distance between the ob- 
servation point and the body. The latter consideration suggests 
that the (iauss-l.egendre quadrature formulation is especially 
well suited for modeling satellite-level gravity and magnetic 
anomalies Iweatise of the large elevations involved with these 
measuremenl,s. The gravity and magnetic anomaly signatures of 
the crustal thickness model, for example, can be computed 
using an iik > it/ x iii = 4 x 4 x 2 - .32 point quadrature formula 
to neaily the same precision as devehiped by (he 512 point 
formula in Fig. 413 and 5, respectively. 
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A|»|riiriiliiNi of \trih«td 

l lic spherical earth iiiodcliiig priicedurc Mas used to ohiain a 
preliminary mcm of the graxily and inapiiciic anomaly charac- 
teristics for the Mississippi embay iiicni at 450km elevation, 
this iiiformaiion is pertineiil to evaliiaiiiip. foi esample. llie 
leasihiliiy of iisiiiy: satellite (ttawly and iii.iyncltv surveys lot 
detecting' aiioiiialy siy>iiattires due to failed tilts. 


The Mississippi emhayment represents a broad, spuun- 
slupcd rc-ciilraiit of Meso/.oic and Ceno7.oie sedimentary rocks 
vvhicti extends into the Paleozoic terrain of the North Amcr* 
iean ctaton from the south as outlined hy the shaded contour 
of Tig.6D. The axis of this feature roughly parallels Ihe Mis- 
sissippi Kiver tapering iiorlh Aard into the tectonically active 
legion of the New Madrid seismic /one. An integrated unutysis 
of gravity, seismic, stratigraphic and petrologic data by Lrvin 
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IS iIk yt.’viiy profii. .iiklis-d h> I rvin aiul Mcliinnis K is the scismiv refraction tine sliHltcd by Msi'amy anil 

surr.icc Have propayalioii palh sluilieil by Atislm aiul Keller (IV7‘t) 


aiiii Me<iir.->iN (1975) suytgcMs the emhaymcni is a late I’re- 
canibriaii aulaeogen which was reactivated most recently in the 
late f’rctaccoux by tcnsioiial forces initiated during the for- 
mation of the prcseiil Atlantic ocean basin by subsidence of llie 
(juir coastal plain. 

t'igure 6.A is a cross-scetion of the density structure of the 
Mississippi cmhayment given by I rvin and MeCiinnis (1975) 
■long a profile between Yellville. Arkan as and Scottsboro. 
Alabama (hereafter called V- S profile). This density model was 
synthesized from regional gravity data derived from the IfS. 
Houguer gravity anomaly map of Woollard and Joesting (I9f>4) 
and the results of a reversed seismic refradion profile between 
l.ittic Kock. .Arkansas and ('ape Girardeau. Missouri as de- 
scribed by Mcl'ainy and Meyer (I9(>6). Austin and Keller 
(1979) integrated the work of MeCamy and Meyer ()966) with 
an analysis of Rayleigh wave dispersion along a propagation 
palh between Oxford. Mississippi and Horissanl. Missouri to 
obtain a similar dciisily model for ihe > - \ prolile which is 
illustrated in l ig. (>1) An index map for locating these various 
studies is given in l-ig.bU. In general, the crustal cross-scctions 
shown in l ig. 6A. B supp<*rt the failed-rift moilcl for the origin 
of the Mississippi embay mciii 

The agreement of surface wave, seismic refraction, and grav- 
ity data in Ihe region of the enibaymeni suggests that the 
crustal cross-section given in Tig. bU can be useful for develop- 
ing ■ reasonably valid three-dimensional misdcl of ihe embuy- 
meni. Accordingly, the crustal cross-scclion that was genera- 
lized from Fig. 6B for the purposes of this study is given as the 
four-body model shown in I ig U I he gravity analysis due to 
Cordell (1977) was used lo proiect tlie characteristics of this 
generalized crustal cross-scetion north and soutli of Ihe ) - .S' 
profile 


Cordell (1977) cv'rreclcd Ihe smoothcHl positive ikoigtier 
anomaly of the cmbaymcnl for the low-density scHlimcnts ami 
observed the long continuous positive anomaly with an ampli- 
tude of 15 45 mgal increasing southward illustrated in l-'ig. fd) 

I he axis of this anomaly closely follows the Mississip|vi River 
norihward beyond its conntiencc with the Ohio River into routhern 
Illinois. The anomaly exhibits relatively uniform behavior south of 
the y--S profile until about 33' N where it increases sharply, thus 
suggesting that the crustal cross-section may be uniformly pro- 
jected southward along the Mississippi River lo approxim.ilcfy 
.1.3'' N. Id the north, decreasing gravity anomaly values in 
conjunction with the northward tapering surface configuration 
of the cnihuyinent suggest a coinmcnsuralc northward tapering 
projcelioii ivf the crustal cross-section along the Mississippi 
River into southern Illinois. Hence, to obtain the first-order. 
Ilircc-ilimensional generuli/iilion of the crustal structure of ihc 
embaymciit u.icd in this investigation. Ihe crustal cross-secli«Mi 
of l ig 6C was projeclcd uniftrrmly south of the l-S profile 
and tapcied uniformly northward as outlined in I'ig. M>. fhe 
nofihern ends of the four bodies of this generalized model as 
projected onto the cross-scetion along the V-S profile are 
given by Ihe shaded regions of F'ig.6C. 

To compute ihc potential field anomalies at 450 km ele- 
vation. each of Ihe four bodies of this generalized model was 
represented by a Gauss- Legendre quadrature formula consist- 
ing of 1 2K equivalent point sources. The latitude and longitude 
limits of each body were represented by X point sources and 
the radial limits by 2 point sources. Pertinent body volume 
limits were iniei|volated from a set of body points that sampled 
the eiHtrdinales of the surface envelope for each bt*dy. Ihe 
qu.idrature formulae were next evaluated and summed over ii 
(2). 1 1) obscivaium gm! spanning Ihe region (2«) 2X0) *1,. (" 
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l-ig. 7A K. Mi'MNMppt cmhiiNiitciii miicIIiU'-IcvcI hciMcco ( I cgcnUrv ()uuJr;iiurc modck'tl ^rsiviiy uml mugnctic ummiuljcti 4nU 

upw.iril coniiiiiicil l°-u\cr;i)!C(t frcc-.iir gravity iitui rvvliivvil i<i-lhv'-pi>lo I'lKiO utlcllilc magttclic ummiiily ( itch iinomiily Ikld it pii'itwi on n 
vlvrcograpliiv t'«|lial-aiv:i |xil;ii |iiok'v'Iioii 


45) ”N ;iiiv) v'oinp.ircil it' oh>cr\ctl gravit) ;ituJ iitagav'liv ;iiu>- 
inaly data at 450 km ck'valion 

t'inv'iiid coiiiimialioii td' Ircv-ait gravily aiiotiial) values 
from the surface of llie carih lo an elevation v*f 450km by 
cv|uivak*nt point source inversion for the study area leads to 
tlic results shown in I ig. 7A Tltese data exhibit a pronounced 
relative positive anomaly with slightly greater than 5 mgal of 
relative amplitude in the region of ihe embay nient the gravity 
effect of the geiierali/ed four-body model described iibove is 
shown in I tg 71) as roughly a 4nig.il anomaly. Ihe general 
agrc*ement K'tween the modeled and observed data over the 
embayment suggests that the observed gravity anomaly can K‘ 
reasonably well accounted fivr at 450km elev.iinm by the gen- 
eralized four-body model. 

IV'l.ir Orbiting Oeophysieal ( Ihseiv.iiory il’(KiO) s.ilellile 
magnetomeier olvservatioirs reduced to i.idial pol.iii/ation using 
.1 norm.ili/.ilion anipltuule v>f (>ll.ilOll g.iiiiiiu by equiv.ileiii |votiii 
source inversion are given for the stiulv area in I ig 7( Ihese 
data show a jvuiminent e.isi-west magnetic high ili.ii is bu-.iched 
m the vicinity of the embayineni by a inagneiiv low lo give 
better resolution of the eh.iracleiistics ol the m.igiieiie .inoinaly 
lor iIk' emb.ivmeiil region, llte ladi.illy polaii/c'd il.n.i woie liigh- 
pass filteied fur anomaly wavelengihs smaller than abotii Id Ihe 
higli-pass filuieil data .iie illiisii ited in I ig 7|> and show ,t neg- 
ative anomaly of roughly .) ganim.i over the enib.iyinenl 

\kasitewski ct al |ld7n| found that most analyses of me 
dium to long-wavelength magnetic anom.dies suggest that sour- 
ces are probably contained in the lower crust which, in gener- 
al. may be stibstaniially more magnetie than the np|vr ertisl 
Ihe conditHnis for coherent regional in.igneii/atioii aie en- 
hanced as crustal depili increases Kcmancnce aiul tlieiinal 


overprints arc diminished, and viscous magnetization and ini- 
tial suseepiihility arc enhanced with increasing tempera- 
ture especially within I00“-I50"(.' of the Curie ptrinl. The 
thickness of the crust within this thermal regime may be 5- 
20 km deiwmdir.g on the steepness of the geothermal gradient. 
Acvordingly, they suggest that deep crustal magnetic sources 
are probably related to lateral variations of petrologic factors 
or Curie isvtiherm topography. 

Acc«>rdingly, an obvktus deep crustal source for the ob- 
served magnetic anomaly is body #2 (I ig. bC) which also re- 
presents the major gravity source of the embayment model. Austin 
and Keller ( I V71) propose that the combination of bodies # I and 
A 2 was formed as a manifestation of a mantle upwurp beneath 
the embayment comprising of a mixture of crust aitd upper 
mantle material which subsc((ucntly cxioied tvv form a block of 
high density material. Magnetie hyptithcses which arc exm- 
sistent with this view include body #2 us ;• zone of negative 
magnetization contrast with respect ti> Ihe lower crust due to 
depletion of magnetic minerals. Negative magncli/ation for 
Utd) # 2 aKo can result from Icinpcralures which exceed the 
Curie |voinl. although present heat How data tSass ti ul 1976) 
do not appear l«« warrant this hypoihc*sis for the embayment. 

UvKiy #.l may rcprc*sent an additional magnetic source for 
the embayment assuming crustal magncli/ation increases with 
depth. However, the pvuilivc magnetic contribution of body 
A .) w ill he relatively weak if the Curie isotherm depth is about 
4Ukm or rnttre. Arguments for including bodies # I and #4 in 
a magnetic model of the embaymeni appear to be lacking, so 
iliai body # 2 probably represents Ihe primary source for the 
observed magiKiic anomaly data if Ihe Curie isotherm is al 
alvout 40 km of depth m Ihe regmn of the embayment. 
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licncc, ihc mii|!nc(K' uniNiiiily due to b«s!) # 3 wii% caUuliil- 
cd ill t ig 71. iMiiig ii magiHMi/iilHm ci>iilni!it of 2.4 
X to 'emu cm*. ihCM.* rcMitiii shim that the unonialy ampli- 
tude iihserveil f»»r the rcgi«»ii of the cmbaymeiit at 450 km 
elevation etui be tiialelied ivell by a vource eueli av KhI> # 2 
tiHutcd near the base of the erunt with magnetic properiicii 
which eor^e^pond with the iiiagiKliration cliiir:tcicrii>;icii anti- 
cipated for the lower ctuM. Subulantiiil disparity ic appurciit, 
howcier. when the spaiiat chariieterixtien of the obM*rved and 
mtHlek'd magnetic anomaliec are comparetl. l-'urther rcniie- 
mcniii of the magnetic model are neeewtary and will be particu- 
larly warrenled when data are available from the current Mag- 
Mit program 1 1 angel IV7tl| to verify and further upgrade the 
IHXtO katellite magnetic anomalies for liiho»pheric appli- 
cations. Accordingly. twHiy ti 2 as determined by gravity 
modeling ci>nsideratioiis represents only a preliminary magnetic 
model for the embay ment. 

C'lNirliMiins 

Regional gravity and magnetic anomaly modeling for arbitrari- 
ly shaped lilhosplierie sources with variable physical properties 
can be achieved accurately and elTiciently m splieric.l c«ho- 
dinates using (ianss-l egendre i|uudruiiire integration. I'he pro- 
cedure involves representation of the anomalous source as a 
distribution of equivalent point gravity pole's or point magnetic 
dipoles which contribute incrementa! aimmaly values evaluated 
and summc'd ut each observation point lo obtain the total 
anomaly. Ihe distribution of c'Muivaleiit point sources is de- 
termined directly from the volume limits ot the anomalous 
body. l-or an arbitrarily shaped body, the variable limits of 
integration can be obtained from interpolations performed on a 
set of body p*nnts which approximate the suiface envelope *»f 
the anomalous s«iurce. 

A chief practical advantage of the methiHl is its consider- 
able versatility. Ihe pliysical pro|ierties of the ev|iiivalent point 
sourevs, for exain|vlc. can Iw individually varied to reflect physi- 
cal property variations of the b«>dy being modeled. The inelh<Hl 
can also be readily extended to model the geoidal anomaly, 
vector components, and spatial derivative's lo any order of the 
biHly's anomalous gravitational and magiK'tic potentials. |-i- 
nally. the accuracy «*f the method can be cmiirolled by adjust- 
ing the number of equivalent point sources or the distance 
between the source and observation point. In this regard, the 
mcihod is particularly well suited for satellite gravity and mug- 
netie Hmmialy modeling K'cause the cfliciency and ac'curaey of 
tin.' application increases with increasing distance betwc'cn 
source and tibscrvaiiou points. 

In consideration of the foregoing results, it is concluded 
that (iauss-l egendre quadrature intcgtaii«m f.icilitatcs a power- 
ful and efficient approach to spherical cariii iiiiidclmg of re- 
gional-scale lilhosplierie gravity and magiK'lic anomaly stuircVs. 
Accordingly, the method has widespread applicaiHin in the 
analvsts and design of rcgional-seate gravity and ni.igiielie sur- 
veys for lilhosplierie Hive's! igalHUi 

I iiiuiisul siipfsirt Cor ihis imcsiigaiHNi was pro- 
vkic'd by the <i»JJaiJ S|Mi.c I liglii ( ciiicr uihlcr K.ASA foniiaci 
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